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Abstract. We develop a perturbation method that generalizes an approach proposed recently to treat 
velocity-dependent quantum-mechanical models. In order to test present approach we apply it to some 
simple trivial and nontrivial examples. 

PACS. 03.65.Ge Solutions of wave equations: bound states 



^ I 1 IntrodllCtion states. Jaghoub chose rather too simple examples in order 
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In a recent paper Jaghoub [Q developed a perturbation 



to test his method £Q. 



^ theory for velocity-dependent quantum-mechanical mod- 



els. Such potentials are useful, for example, for the study 
of pion-nucleon scattering and in models of particles with 
coordinate-dependent masses [1] (and references therein). 
According to the author, one of the advantages of the 
method is that its main equations depend only on the 
eigenfunction of the chosen unperturbed state. More pre- Th e purpose of this paper is to generalize the method 



Similar perturbation approaches are known since long 
ago, but they apply mainly to local (coordinate-dependent) 
potentials [2]. Although one can in principle adapt some 
of the well-known perturbation methods 3 to the treat- 
ment of velocity-dependent potentials, here we proceed in 
a different way. 



cisely, the method does not need neither the whole unper- proposed by Jaghoub Q] and derive a perturbation algo- 
turbed energy spectrum nor the basis set of unperturbed rithm for the treatment of a wider variety of problems. In 
eigenfunctions as in the formulation of the Rayleigh-Schr6din|§?-Elwe develop our version of the method in a quite gen- 
perturbation theory in terms of sums over intermediate eral way. In Sec.|3|we apply the perturbation approach to 



Correspondence to: F. M. F. an exactly solvable example which enables us to comparc 
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the results of present procedure with the expansion of the 
exact eigenfunctions an eigenvalues. In Sec.^we discuss a 
partially solvable example treated by Jaghoub [T . Finally, 
in Sec.Elwe summarize our main results and discuss other 
potential applications of present method. 

2 Method 

Suppose that we want to solve the differential equation 



y"(x) = F(y,x) 



d) 



where the prime denotes differentiation with respect to x 
and F(y,x) an arbitrary linear differential operation on 
y{x). In order to apply perturbation theory we choose a 
closely related, solvable problem 



Vo( x ) = F o{yo,x). 



(2) 



If we multiply Eq. Q) by yo and subtract Eq. (J2J) multi- 
plied by y we obtain 
d 



dx 



W(y, y ) = Fy - F y 



where W(y, yo) stands for the Wronskian 



W(y, y ) = y'yo - y' y = y 



2 ( y_ 

yo 



(3) 



(4) 



On integrating Eq. twice we obtain 

f x dx' 

y{x) = C 2 y {x) + C 1 y {x) / - + 

J p yo{x'y 

J 13 yo[x) Ja 

where the integration constants C\ and Ci and the in- 
tegration limits a and /3 enable one to accommodate the 
boundary conditions of the problem and the normalization 
of the solution. 



Notice that the function 



u(x) = y (x) 



dx' 



yo(x') 2 



(G) 



satisfies W(u,yo) = 1, so that u(x) does not vanish and is 
finite at the zeroes of yo (x) . If the operation Fo is simply of 
the form Fo(yo,x) — fo(x)yo{x), then yo{x) and u{x) are 
two linearly independent solutions of Eq. J2J), and u{x) is 
called ghost state |1| (and references therein). A complex 
linear combination of yo{x) and u(x) proved suitable for 
the construction of a logarithmic perturbation method for 
excited states 0]. 

In order to apply perturbation theory we introduce a 
perturbation parameter A into F and expand 



and 



F = J2FA j 

3=0 



3=0 



(X) 



(8) 



We simply introduce the series l|7| and (JSJ) into Eq. (JSJ 
and obtain y j (x) in terms of yk(x), k = 0, 1, . . . , j — 1, and 
Ek, k — 1,2,..., j. The boundary conditions determine 
the appropriate value of E j that is the only unknown in 
the expression of y j (x) . The procedure will be made more 
explicit in the examples below. 

Present approach applies to bound and unbound states; 
in this paper we concentrate on the former ones. It is cus- 
tomary to choose a convenient normalization for bound 
states. Here we keep our equations as simple as possible 
and add a normalization factor N(X) = Nq + Ni\ + . . . 
at the end of the calculation if necessary. We may, for ex- 
ample, require that (A^ + A^A + . . .)(yo + J/iA + . . .) be 
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normalized to unity up to a given perturbation order. As- being the exact form of the normalization factor A is un- 
suming that yo(x) is normalized to unity we easily prove necessary for present perturbation calculation and we can 
that add it easily at the end of the process as discussed above. 

In order to simplify the presentation and discussion of the 



A 2 

N(X) = 1 -A / y yidx + — 



VoVidx ) (9) 



results we arbitrarily choose the normalization factor for 
-2 J y y 2 dx - J y\dx +..., (10) the unpe rturbed solution: A = A(X = 0) = y/2. 

Straightforward expansion of the exact results 112(1 in 
a Taylor series about A = yields 



where the integrals extend over the whole physical coor- 
dinate range. Obviously, the eigenvalues are independcnt 



x- 



of the normalization constant. yj{x) — . yo(x) 

It is worth noticing that we have not made explicit E (n) = n 2 TT 2 5 o + — S 2 (13) 

use of the assumed linearity of F(y, x) in order to derive 
Eq. (JSJ). Consequently, that equation applies to arbitrary 
nonlinear differential operators F and Fq. However, the 
eigenvalues of nonlinear equations do depend on the cho- 
sen normalization condition and we have to take this fact 
explicitly into account in order to solve them. 

3 Solvable example 

In order to test the general equations developed above we 
choose the eigenvalue problem 



In order to apply present perturbation theory to this sim- 
ple test model we make the obvious choice F = Xy'(x) — 
Ey{x) and Fq = —Eoyo(x), where Eo(n) = n 2 ir 2 and 
yo(n,x) = v / 2sin(n7rx). Since u(x) does not vanish at the 
end points we choose C\ = 0. The other integration con- 
stant C2 is important to force a normalization condition 
at each order of the perturbation algorithm. However, as 
in this first step we are not interested in the normalization 
of the solution we simply set C2 = 0. Thus we are left with 
the hierarchical perturbation equations 



f x dx' f x 

yj(x)=yo(x) / (yoy'j-i (14) 

y"(x) = Xy'{x)-Ey(x), ^ yo(x') 2 J 



V(0)=y(l)=0 (11) 

with solutions 

y{x) = Ae Xx/2 sin(n7ra;) 
,22 A 



-y J2Ekyj-k)(x")dx", (15) 
fe=i / 

where we have chosen a — [3 — to satisfy the boundary 
condition at x = 0. Notice that yj depends on Ej that 
is determined by the boundary condition at a; = 1. For 



E — n 7r + — , n — 1, 2, . . . (12) example, at first order we obtain 

x cos(ri7ra; ) / \/2x x/2E 1 \ , . 

Notice that the eigenvalue equation 111(1 is not Hcrmitian yi(x) = h sm(n7rx) 

2n7r \ 2 2n A Tr 2 I 

but it supports real eigenvalues for real A. For the time (16) 



4 Paolo Amore, Francisco M. Fernandez: Wronskian perturbation theory 

that does not satisfy the boundary condition at x = 1 mi- / _ nirx(x 2 — 1) cos(n7rx) / 3x 2 1 . 

y\\x) — — h i — — — h -rr l svQ.\nirx), 



10 V 20 10 

less E 1 ! = 0. Proceeding exactly in the same way we obtain 3(8n 4 7r 4 + 10n 2 n 2 — 15) 



E- 



2 



1000n 2 7r 2 



the function and energy perturbation coefficients (1131) or- , , 

5J r 248n 6 7r 6 + 462n 4 7r 4 - 1575n 2 7r 2 + 1890 , on . 

der by order, which clearly shows that Eq. (|15(l gives the 35000n tt 

correct answer. The energy coefficients agree with those shown numeri- 

The normalization factor calculated by perturbation ca n y by Jaghoub 1 for n = 1. In this case the normaliza- 

theory tion factor is given by 



_ , A A 2 (n 2 7r 2 + 12) A 3 (n 2 7r 2 -12) 

iV(A) - 1 "4 + 96n^ 2 + 384n 2 7r 2 + ' ' ' ^(A) = V2 



1 3A 3(29n 4 7r 4 + 10n 2 7r 2 + 30)A 2 



also agrees with the one derived from the exact solution. 

4 Partially solvable example 

In order to compare our approach with Jaghoub's one 
more closely, in what follows we treat one of that author's 
examples: 



20 4000n 4 7r 4 

(21) 

We do not show the analytical expressions of 2/2(2) and 
2/3(0:), as well as perturbation corrections of larger order, 
because they are rather complicated. 

5 Conclusions 



3x 2 \ ... 6 r ... ... „ . N _ In this paper we have shown how to develop a useful per- 

1 - — y"{x) - - [xy'(x) - y{x)} + Ey{x) = 0, (18) 



5 / 5 

turbation method from the Wronskian between the per- 
where y(0) = 2/(1) = exactly as in the preceding exam- 

turbed and unperturbed solutions. The approach enables 

ple. Jaghoub constructed this problem in order to have an 

one to treat velocity-dependent quantum-mechanical prob- 

exact solution for the ground state 

lems as well as local perturbation potentials. Present ap- 
y(x) = Ax(l — x 2 ), proximation is more general than one proposed earlier pQ 

E(n — 1) — 6 (19) anc ^ ename s us to diseuss mathematical aspects of the so- 

lutions that could be otherwise masked. We easily derive 
the most general equation for the application of perturba- 
tion theory and can clearly analyze all the contributions to 
the approximate solutions and their behaviour regarding 
normalization and boundary conditions. 

We have tested our general equations on an exactly 
solvable model showing that the approximate method yields 

2n 2 7r 2 + 15 

Ei — — , exactly the same results that one obtains from expansion 



and treated only this trivial case by perturbation the- 
ory. Here, on the other hand, we also consider the ex- 
cited the states for which there are no exact solutions 
as far as we know. A straightforward caleulation with 
2/0(2) = sin(n7ra), F = 3x 2 y"/5 + 6(xy' — y) /5 — Ey, 
and .Fo = — £-o2/o shows that 
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of the exact solutions. We have also treated a partially 
solvable problem. In both cases we have derived perturba- 
tion corrections for all the states in terms of the quantum 
number, thus generalizing Jaghoub's results for the latter 
model p. 

Present perturbation method also applies to nonlin- 
ear models. We may treat them exactly as the examples 
above, except that one has to consider the normalization 
condition explicitly at every perturbation order because 
the eigenvalues depend on it. 

For simplicity we have chosen the unperturbed and 
perturbed equations in such a way that we could solve 
all the integrals analytically. In more difhcult cases one 
should have to resort to numerical integration. However, 
any problem that can be treated exactly by Jaghoub's ap- 
proach can also be treated exactly by present one, proba- 
bly in a more general way. 

By means of an appropriate choice of the boundary 
conditions present perturbation approach is suitable for 
the approximate calculation of scattering phase shifts. In 
this field our method may be an alternative to the loga- 
rithmic perturbation theory that is commonly applied to 
local potentials 
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